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I. INTRODUCTION 

The metric for a stationary cylindrically symmetric spacetime has the Papapetrou form [Q 

ds 2 = -e 2U (dt + Ad V f + e- 2U [e 2K (e 2H dr 2 + dz 2 ) + W 2 d V 2 } , (1) 

where the metric functions depend only on the radial coordinate r. The function H determines the gauge and can 
always be scaled to the isotropic form H = by changing r. The static case follows when A = 0. When the source of 
the gravitational field is a perfect fluid, the energy-momentum tensor reads 

T a/3 = (fJ, + p) u a up + pg af 3, (2) 

where p, is the energy density, p is the pressure and u a is the four-velocity of the fluid. Rigid rotation is characterized 
by vanishing shear. The field equations are less than the unknowns, which allows to prescribe an equation of state 
p = p(p). 

The form (2) vanishes in vacuum, p = fj, = 0. Dust is another particular case with p ~ 0. Vacuum solutions with 
a cosmological constant A are obtained when p = — p = A/k. Here k is the Einstein constant. When p = — A/k, 
p = p + A/k we have dust solutions with density p on a cosmological background. 

The first solutions with a source given by Eq. (2) were found by Lanczos in 1924 ||. He solved the case of dust 
and discussed the case of dust plus a cosmological constant. Several years later the general vacuum solution was 
found by Lewis g]. The first global solution, containing rotating dust interior, was given by van Stockum |j who 
rediscovered the dust solution of Lanczos. A rotating universe with cosmological constant and constant p = p = —A/k 
was obtained by Godel ||. A would be new five-parameter solution || was shown much later j?J to be equivalent to 
the two-parameter Lanczos solution. 

The field equations for perfect fluids with non-zero pressure in the axially symmetric case were first written by 
Triimper ||]. The cylindrical problem was studied by Krasinski who obtained six types of new solutions, all but one 
expressed in terms of hypergeometric functions J9|-[l^] ■ The equations of state of these models are rather complicated. 
Meanwhile, in the static case a solution with the linear equation of state 

P={l- l )^+Pi, (3) 

where 7 = 6/5, was found and the general problem with p\ = was simplified Jl3|. The case 7 = 4/3, pi = was 
solved next |l4| . The general solution with vanishing p\ was soon derived |l5| and rediscovered in a more convenient 
form some time later [^6|. The connection between the two was clarified in Ref. fTij] . The static solution with 
cosmological term, which corresponds to 7 = 0, p\ — was derived too p8[ . General static solutions with p\ 7^ 0, 
however, are accessible even today only numerically, including the simplest case of constant /1 [ fl9| . 

The stationary metrics were studied in the past 10 years along two basic lines. The general cosmological solution, 
which is implicitly contained in the Krasinski solutions |2(| , was obtained explicitly |2l]-|2j| . This made it possible to 
build an exterior for the Godel metric p4| . On the other hand, solutions with pi = and 7 in the interval allowed 
by the energy conditions (1,2] were found by trial and error by Davidson [25|-|2^|. They confirm the results of a 
qualitative numerical study, using a dynamical systems approach pSj]. Davidson discovered by the same method a 
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solution with finite radius having 7 = 2/3 and p\ ^ p{| which he matched to the external Lewis solution [ j30[ . 
Finally, Sklavenites studied two models with complicated equations of state, which reduce to the Kramer and Evans 
solutions in the static limit Q. 

In this paper we find analytically the gravitational field of a rigidly rotating perfect fluid cylinder with equation of 
state given by Eq. (3) with p\ = 0. Connections are made with most of the results found in the past. 

In Sec. II the system of field equations from Ref. jlj is simplified and the regularity conditions are discussed. In 
Sec. Ill it is solved analytically after a proper choice for H . The interval 1 < 7 < 2 is divided into four regions. Two 
of them possess realistic algebraic solutions. They are separated by the exceptional point 7 = 3/2 where the solution 
becomes exponential. 

In Sec. IV some particular solutions, mentioned in the Introduction, are obtained from the general solution when p 
and fj, are constant. The dust solutions in vacuum and in a cosmological background are rederived for completeness. 

The properties of the general solution are studied in Sec.V. The acceleration, the vorticity vector and the angular 
velocity of rotation of the perfect fluid are found. The curvature invariants are regular. Closed timelike curves (CTC) 
appear in all cases. Sec. VI contains some conclusions. 

II. FIELD EQUATIONS AND REGULARITY CONDITIONS 

We shall use the system in Ref. jjj which is in comoving coordinates. Thus, u* = e~ u , the other components vanish. 
All dependence on H may be hidden in the operator ' = e~ H -4-, which possesses many of the basic properties of a 
derivative. Then the field equations read 

A'W 

A" + AA'U' - = 0, (4) 

= ltope* K - u \ (5) 

„4U 4/2 TT'W' 

U" + e -^ + ^=4.(» + 3 P )e^-«\ (6) 
W" ,W ,0 e AU A' 2 

K" + U' 2 + ^^=8irpe 2 ( K - u \ (8) 

p' + (fj, + p)U' = 0. (9) 

We use units where G = c = 1. Then k = Sn. The last equation follows from the others. Together with the equation 
of state (3) this is a system of 6 independent differential equations for 7 unknowns \x,p 1 U, K, A, W, H . The solution 
should depend on the arbitrary gauge function H . A regular solution must be elementary flat. This means 

lim r- 1 e u - K - H ( e - 2U W 2 - e 2U A 2 ) 1 ' 2 = 1. (10) 

Let us simplify the above system. Eq.(4) is easily integrated once 

A' = a We' 4U , (11) 

where ai are constants. This holds for any equation of state and simplifies the terms with A in the other equations. 
When 7 f 1 Eqs. (3) and (9) give 

p = p e ^ , (i= -e T-i . (12) 

7- 1 

For stability reasons the constant po > 0. It measures the central pressure. Eq. (6) can be integrated next 
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(13) 



Eqs. (7,8) give an equation for K 



WK' -W = -yA + a 2 . (14) 



Eqs. (13,14) show that K — U may be expressed through W. We have reduced the initial system to Eqs. (5,11-14), 
all but one being of first order. 

In a realistic solution U, K and H should approach constants when r — > 0. Then Eq. (10) shows that W <~ r v ,A ~ 
r T . However, Eq. (11) gives r = v + 1. Therefore, the term with A decouples in this limit and we get the static case 
condition 



Hence, W ~ r and A may be written as 



where q is a constant and F (0) = 0. 



l ime -*-"M = i. (is) 

j— >0 r 



(16) 



1 + F (r) ' 



III. THE GENERAL SOLUTION 

The key to the solution of the simplified system of equations is the choice of gauge. Different gauges have been 
proposed in the literature but none of them works well here. Let us choose the following H 

h _ aohW r 

16ir P o(A + ho)' 

where h, h are constants. The introduction of ho is necessary to ensure the regularity of e H on the axis of rotation. 
We also require H (0) = 0. It is possible to satisfy Eq.(17) by fixing appropriately most of the constants without 
fixing H . Therefore the solution will depend on an arbitrary function. 
Eq.(5) becomes a relation between K and U 

he 2K = e^ u , (18) 

which shows that h > 0. Then Eqs. (13,14) should be checked for compatibility. This is achieved by adjusting ai and 
h 

2 - 7 4 - 37 , . 

° 2 = 2(7" Yf 1 + I^t - lf oV (19) 

. 8tt Po II7 2 - 247+12 

^^r /(7) ' /(7)= (7-D07-4) • (20) 

The zeroes of the nominator of / are 71,2 = 12±|y2_ xhe second is outside the realistic range, while the first is 

12 + 2^ 2V3 



11 2-s/3 — 1 



1.4065. (21) 



The sign of / divides the range of 7 into three regions: (1, 4/3) and (1.41, 2) where / > and (4/3, 1.41) where / < 0. 
The sign of po must coincide with the sign of /, therefore in the third region the pressure is negative and unrealistic. 
Also / blows up at 7 = 4/3. We can arrange for U and H to vanish on the axis. Eq.(18) then shows that the same is 
possible for K only when h = 1. This condition is a relation between the two independent parameters cto and po of 
the solution 
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"f-S- < 22 > 

It also means that K ~ [/ '. These restrictions characterize the Davidson solutions. 

In general, e x — * h~ x / 2 and — > hr x l 2 r (from Eq. (15)) on the axis. Then Eq. (17) fixes ho while Eq. (11) yields 

1 



Let us introduce 



1/2/ sz.gna a , . 

8 V 7T p s«5 n Po 2V/i 



F= (A + h ) 2 . (24) 



Inserting Eq. (17) into Eq. (11) and taking W as an independent variable, we get 

■*~^ + ir < 25 > 

The last equation to be satisfied is (13). Taking into account Eqs. (17,25) it becomes a second-order non-linear 
ordinary differential equation for Y (W): 

, - ; , C \Y W 

I WW 



(^ + tW- < 26 > 



i = |l^, Co = - / (H£l + M . ,27) 



37 — 4 \ ao 

The parameter Co determines 01. Eq. (26) is a generalized Emden- Fowler equation [^T). It belongs to the class 2.6.2.98 
and, fortunately, is among the 147 integrable cases. The parametric solution hints at the possibility to reverse the 
variables. Writing W = cie G , where c\ is another constant, one obtains an integrable equation for G which yields 

G = I c 2 + (b + l)Y + 2c Y^ (28) 
with yet another free constant ci- This determines W as a function of A. Eq. (17) becomes 



6 =(2ao) ~ATho~' (29) 



The parameter ao measures the effects of the fluid's rotation through Eq. (16) 

A= ■ (a ° r ? 2 r. (30) 

2v^o7(l + F) v ; 

The static case is given by ao = 0. Obviously Eq. (17) or (29) is tied to the stationary case. In principle, A 
may be expressed through H but an integration is necessary, which depends on F. It is much easier to give all fluid 
characteristics as functions of A. Next we determine U from Eqs. (11), (29) 

fWW A 

2{A+h Q y {31) 

The function K follows from Eq. (18) and is linear in U . The free parameters for each 7 are two: ao (measures 
rotation) and po (measures the central pressure and density /io = po/ (7 — 1)). The Davidson solutions have just one 
parameter because po and ao are related by Eq.(22) like the mass and the angular momentum of the Kerr solution. 

The integral in Eq. (28) has two cases: 60 + 1 = (7 = 3/2) and bo + 1 7^ (7 7^ 3/2). In the second case it is given 
by two expressions, depending on the sign of Cg — (60 + 1) C2- The condition W ~ r near the axis is very restrictive 
and tedious calculations show that W looses its dependence on cq and C2 
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-i /2 
W = W A 1 / 2 B*&=*!, |W |= / 1/a > (32) 



S=l + aA, a= ,f 7 3 > ■ (33) 
(37 - 4) ft 



We can always choose a > 0, q > and Wo > so that A > 0, B > 0. The metric functions and the pressure read 

VP 



e H = — 2(2^-3) ) e 2 ' 7 = S (34) 



7(2-7) „„ , (2— ,y 

p = p a B *(t-i)(27-3) , e 2A = /i" 1 B 4 (^~ 1 )< 2 ^- 3 ) . (35) 
It is useful to write the metric in the so-called Lewis form, used in most of the research on the topic 

9 2-t 2-7 2- T j4? 97 2 -227+12 

ds 2 = -B^-^dr - 2AB 2 < 2 ->- : » dtdtp + (WqB - A) AB 2 < 2 ^- 3 ) d^ 2 + ^f-j B 4 < 2 ^~ 3 >^" 1 > dr 2 + 

(2-7)(37-4) 

+/l" 1 J B 4(27-3)(7-D dz 2 . (36) 

In the process of solution the gauge function was shifted effectively from H to A and then to F. The simplest 
function should be chosen. However, it is not possible to do this simultaneously in the whole region 1 < 7 < 2. It is 
divided into four parts. 

When 3/2 < 7 < 2 we can choose F — 0. Now / > 0, a > 0,pa > 0, hence B — > 00, p — > at infinity and that is 
physically realistic. Setting h = 1 we obtain the solution from Ref. |^6| or case a) from Ref. |27| ]. 

When 1.41 < 7 < 3/2, / > 0,po > but a < and if F = 0, B will vanish at some point, resulting in a singular 

metric. Let us take F = —ar 2 , then B = (l — ar 2 ) 1 is always positive and decreases monotonically as r — * 00. Then 
p — > as necessary. The subcase ft, = 1 leads to case c) from Ref. (27) . 

In the region 4/3 < 7 < 1.41, / < and to ensure the positivity of h the central pressure should be negative. Then 
Eqs. (17), (30) make sense. The appearance of tension, however, leads to instability with respect to perturbations. 
When h = 1 this case was discussed in |27j . It seems that no physical solution exists for this interval. 

Finally, in the region 1 < 7 < 4/3, f,Po,a and A/h are positive. Then B > 1 according to Eq. (33). B should 
decrease monotonically with the distance in order that p also decreases, but this is not possible since B (0) = 1. The 
solution can not have a regular axis, in accord with the results of Refs. ^j,^8|. One can find a solution whose other 
features are regular by setting A ~ (l + /c 2 r 2 ) 1 as was done in case d) from Ref. [ p7| . 

The general conclusion is that a regular solution exists only in the relativistic interval 1.41 < 7 < 2 and two different 
choices of gauge are necessary to cover it on both sides of the exceptional point 7 = 3/2. The Davidson one-parameter 
metrics are obtained when one puts the correct F and h = 1 into Eq. (36). 

In the special case 7 = 3/2 the integral in Eq. (28) is simpler and elementary flatness requires 



W = b 1 A^e" A , a=J 8 -^, h = ^-, (37) 

V 3 a^fq 



where b\ is determined, like Wq before, by the condition H (0) = 0. Then 



e H = J^ e <*A_ / 38) 



Eq (31) determines U and p, (i, K follow from it 



e 2U = ^V3^bxe aA , (39) 



e 2K^^ e u^ p=poe -X' i ^ = 2 p. (40) 
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The basic difference from the general case is the change B — > e aA . The simplest choice for A = qr 2 leads to a regular 
solution with decreasing pressure towards infinity. It also possesses two parameters. Davidson proposed a single 
parameter solution with A — qr and h = 1 M] , but its gauge function blows up at the axis and it is not elementary 
flat. One obtains 2 on the R.H.S. of Eq. (15). 

It should be stressed that the static subcase is not a limiting case of the general solution. Eqs. (4,11) do not hold 
and Eq. (17) is meaningless. One should start with the system (5-9) with A = and solve it a new. Bronnikov 
used the gauge e H — W, while Kramer p6[ expressed the metric components as functions of U, which can be derived 
from Bronnikov's solution by changing the gauge [p"7f . 

IV. SOME PARTICULAR SOLUTIONS 

The previous formulas are inapplicable when 7 = 1 and the pressure vanishes. Then Eq. (5) shows that one can 
choose W = r. No exotic gauge is needed, so that H = 0. Eq. (9) becomes fill' = and comprises two cases. When 
fi = we have vacuum and the system of Eqs. (4-8,11) gives 

A' = a re- 4U , — = 2U' 2 - fe~ w , (41) 
r 2 

U" + — + ^e- 4U = 0. (42) 
r 2 

Eq. (42) coincides with Eq. (A4) from Ref . |32| , one of the few papers where the Papapetrou form of the metric is 
used. 

In the dust case /j, ^ and consequently e u = 1. Then 

A=|r 2 , K = -fr 2 , M=f|e^, (43) 

which is the Lanczos solution |^,|| . 

Let us discuss next solutions with constant p ^ and [i. If Eq. (3) holds then Eq. (35) holds too and p, /i are 
constant when 7 = or 7 = 2, i.e. p = ±/i. On the contrary, when p, [i are constant (fi + p)U' = so that either 
p = —[i (cosmological solutions) or U = and Eq.(6) shows that K = const = 0. Then Eq.(ll) simplifies. The 
comparison of Eqs. (5,7) yields a 2 = i2-np. When inserted in Eq. (6) this gives p = \x. This is the Godel solution. 
Let us derive these two solutions from the general formulas (32-36). The cosmological constant is A = — 87rpo in all 
cases that follow. 

a) Godel solution. We have 7 = 2, p = po = u, f = 4, U = K = and 



32wp a 



B = 1 + aA, a= ^2np 0: (45) 

W = WqVAB, Wl = ho- (46) 
The gauge is again H = 0. This implies an equation for A: 



A r = 2y/qAB. (47) 



Its solution reads 



A 1 ' 2 = 1 ( e lr+c - a- x e- lr - c ) , (48) 

where I = ao/2%/2, while c is an integration constant. It is fixed by the condition A (0) =0 to e c = l/y/a and the 
solution has two parameters ao and po. It is an elementary flat subcase of the Wright solution ||, Eq.(31). Then 
hyperbolic functions are obtained. If we fix the two parameters in terms of one constant b 
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the one-parameter Godel solution follows |2J] 

ds 2 = -dt 2 - 4^26 sinh 2 ^rdtdip + 4b 2 (sinh 2 ^- - sinh 4 dip 2 + dr 2 + dz 2 . (50) 

b) Cosmological solution. We have 7 = 0, p = po = —fi, / = 3 and 



% V 327rpo 



B = 1 + aA, a = ^/ 67rpo , I = y aqh — a, (52) 
W = WoA 1 / 2 ^ 1 / 6 , Wo 2 = -^=, (53) 

VD7TP0 

e 2U = B- X l\ e 2K = h-'B 1 ' 3 . (54) 
Let us fix the gauge demanding H = U — K . The equation for A is very similar to Eq. (47) 



A r = 2 v / qhAB. (55) 

Its solution is given by Eq. (48) with a, I taken from Eq. (52), c is fixed as above. The solution has two parameters but 
oq enters only in e 2K . Comparison can be made with Refs. ]2l|-|23[] after the calculation of the quantity ^ = We K ~ u 



* = ^-f = d cosh Ir + C 2 sinh Ir. (56) 

Here Ci are combinations of po,ao- The second relation coincides with Eq. (38) from Ref. ^T|. The case of A with 
opposite sign is obtained by taking negative p$. This makes / purely imaginary and the hyperbolic functions become 
trigonometric. 

Finally let us discuss dust solutions on a cosmological background. The pressure is constant while fi is determined 
from the equations. We set again H = 0. Eq. (9) gives U — and then Eq. (11) becomes 

A' = a W. (57) 

Replaced into Eq. (6) it turns into an expression for fi as a function of K 

^ = -3 Po + ^e- 2K . (58) 
Eq. (8) becomes an equation for K which may be integrated once 

K' 2 = 87r Po e 2K -^K + 2d 1 . (59) 

Here d\ is an integration constant. The term with po may be expressed through K and W in another way from Eqs. 
(5,7). The comparison with the last equation yields 

W = d 2 K', A = a d 2 K + d 3 . (60) 

Everything is expressed through K that satisfies the non- linear first order differential equation (59). It was first found 
by Lanczos and can not be solved analytically. However, we can change the variable from r to y where K = — y 2 . 
Then 

2K 2 4e~ 2 y 2 y 2 dy 2 
e dr = — 2 (61) 
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and all metric components become functions of y |q] . Elementary flatness fixes three of the constants 

V2 

di = -47rpo, d 2 = — , d 3 =0 (62) 
a 

so that the solution has two parameters po.ao 0, like the general solution and the other particular solutions discussed 
in this section. 

This problem was discussed again recently ]33]] but the system of equations was not transcribed correctly from Ref. 
|p4| and the wrong conclusion that an explicit solution is impossible was stated. The system in Iftime's Eq. (17) is 
completely different from Eqs. (4-9) or Eqs. (57-60), the main error being in the analog of Eq. (57). The system 
(4-9) is equivalent to that of Trumper || and its complex reformulation Jm) . We have checked this system also by 
GRTensor. 

In conclusion, we have shown how many of the known particular solutions follow from the general one or directly 
from the system (4-9). Fluids with p\ in Eq.(3) have finite radius. The main system may be simplified again, but 
we have been unable to solve the master differential equation. 

The properties of these special metrics have been studied extensively through the years. Therefore, in the next 
section only the properties of the general solution will be discussed. 



V. PROPERTIES OF THE GENERAL SOLUTION 



We have calculated with the help of GRTensor the characteristics of u a in the metric (1). The expansion and shear 
vanish as it should be for a rigid rotation. The acceleration has only a radial component v r and its magnitude is 

v=(e u ) r e- K - H . (63) 

The vorticity vector has only a longitudinal component lo z and a magnitude 

. = ^e— (64) 
The angular velocity = uj (0) of the rigid-body rotation of the fluid is 



fl = qh = y/2npof, (65) 

where Eqs. (15-18), (20) were used. Amazingly, ao, which determines the magnitude of A and the degree of station- 
arity, does not enter this formula. The angular velocity is entirely determined by the central pressure. This issue 
is obscured in the Davidson's solutions because of the relation (22). Of course, ao has its influence upon the radial 
distribution of v and lo. 

For the general solution these characteristics are 
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3 7 -4 



2npA 



V^Pf- (66) 



When h = 1 and Eq. (35) is taken into account we obtain the expressions of Davidson, his Eqs. (5.7) and (5.9) 
from [^6| or Eqs. (2.16), (2.17) from f27j]. The vector v a vanishes both at the axis and at infinity when F is chosen 
accordingly and 1.41 < 7 < 2. At the axis u> equals fi while at infinity it vanishes together with the pressure. The 
spin of the fluid vanishes at 7 = 1.41 because this is a zero of / (7). 
In the special case 7 = 3/2 Eq. (65) still holds while 

v = kiA l ' 2 e-l aA , w = k 2 e~i aA 7 (67) 

where fcj are constants depending on ao,po- The behaviour of v and ui is as in the general case. 

We have studied also the behaviour of the independent curvature invariants. They are regular at the axis. For a 
typical quadratic invariant the most dangerous term at infinity has the form 
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Fortunately, / ~ v 2 and vanishes at infinity together with v. The linear invariant R has terms I ~ v and so on. The 
Petrov type of the general solution is I except for 7 = 2 where it is D J2(|^7j . 

Finally, let us discuss the appearance of CTC. This happens when g w becomes negative. From Eq. (1) and the 
expression for the general solution it follows that 

g w = e- 2U (W 2 - e AU A 2 ) = AB^J (w 2 . (69) 

The term in the last bracket is always negative for large distances. When 1.5 < 7 < 2 we set F — and g vv < if 

r 2 > 4(7 ~ 1}/ . (70) 



7«5 

When 1.41 < 7 < 1.5 we set F = -ar 2 and g vv < if 



r>— . (71) 



When 7 = 1.5 the second equality in Eq. (69) does not hold. We have 

W-e^A^b.A^e^ (l-^A^j , (72) 

which turns negative if 



A>^, r> 2 -^. (73) 



Thus, CTC appear in all cases and the bigger ao, the closer to the axis they stand. 



VI. CONCLUSIONS 



In this paper we have found analytically the gravitational field of a rigidly rotating perfect fluid cylinder with 
equation of state given by Eq. (3) with p\ = 0. The system of differential equations following from the Papapetrou 
form of the metric is simpler than the system corresponding to the Lewis form. The general stationary solution is 
similar to the Davidson solutions and to the general static solution. It is easier to challenge rotation than the case 
Pi 7^ in the equation of state. The solution is algebraic except at 7 = 3/2 where it becomes exponential. It has 
two free parameters po and ao describing the central pressure and rotational effects. The Davidson one-parameter 
solutions are obtained when Eq. (22) is imposed on the parameters. 

The crucial technical point is the choice of H, valid only in the stationary case. This is a partial fixing of the gauge, 
because an arbitrary function still remains in the solution. It is advantageous to shift the gauge function from H to A 
and F. The two realistic regions of 7 are (1.41, 1.5) and [1.5, 2]. A regular solution requires different choices of F in 
them. Therefore, the gauge should be fixed completely as a last step. As a whole, our results confirm the qualitative 
picture drawn in Refs. p5|-p8|. It is strange that realistic solutions exist only for highly relativistic perfect fluid and 
there is a gap between it and dust. 

We have shown that the cosmological and the Godel solutions follow directly from the general solution. The dust 
solution is easily restored from the original system of equations. The same is true for dust on cosmological background. 
We have tried to solve the controversy between a recent paper and the older research on this topic. 

Interesting features of the general solution are the following. The magnitude of the vorticity vector is a square root 
of the pressure. As a consequence, the angular velocity depends only on the central pressure and not on ao. The 
curvature invariants are well-behaved because the most dangerous terms are proportional to powers of the acceleration 
vector's magnitude. The latter decrease monotonically and vanishes at infinity. CTC appear for any 7 inside the 
realistic region. The parameter ao commands the distance at which g vv becomes negative. More vigorous rotation 
invokes CTC nearer to the axis. 
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